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ABSTRACT 

Motivated by the relative lack of neutral hydrogen around Lyman Break Galaxies 
deduced from recent observations, we investigate the properties of the Lya forest 
around high redshift galaxies. The study is based on improved numerical SPH simu- 
lations implementing, in addition to standard processes, a new scheme for multiphase 
and outflow physics description. Although on large scales our simulations reproduce 
a number of statistical properties of the IGM (because of the small filling factor of 
shock-heated gas), they underpredict the Lya optical depth decrease inside 1 Mpc h^^ 
of the galaxies by a factor of « 2. We interpret this result as due to the combined 
effect of infall occurring along the filaments, which prevents efficient halo gas clearing 
by the outflow, and the insufficient increase of (coUisional) hydrogen ionization pro- 
duced by the temperature increase inside the hot, outflow-carved bubble. Unless an 
observational selection bias is present, we speculate that local photoionization could 
be the only viable explanation to solve the puzzle. 
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1 INTRODUCTION 

Galaxies form from the intergalactic medium (IGM), pro- 
cess such gas into stars, and possibly re-eject a fraction 
of it, enriched by nucleosynthetic products, back into the 
intergalactic space via powerful supernova- driven outflows 
(Mac Low & Ferrara 1999; Ferrara, Pettini & Shchekinov 
2000; Madau, Ferrara & Rees 2001; Scannapieco, Ferrara 
& Madau 2002; Theuns et al. 2002). The energy deposition 
connected to these processes is expected to leave at least 
some detectable imprints on the physical state of the IGM. 
Thus, it is conceivable that such signatures can be stud- 
ied through QSO absorption line experiments. Naively, the 
presence of hot outflowing gas should result primarily in two 
effects: [i] a decrease of the gas density and [ii] an increase 
of the temperature caused by shock-heating (acting in con- 
junction with photo-heating by the UV background) in a 
large (several hundreds kpc) region around the perturbing 
galaxy. Both these occurrences would imply an increasingly 
more transparent Lya forest when approaching the galaxy, 
i.e. a galactic proximity effect. Quantitative confirmation of 
this scenario has faced tremendous difficulties, standing the 
complications of the physics of star formation, explosions 
and metal mixing in multiphase media. Hence, most simula- 



tions to date had to rely on ad hoc recipes for such processes. 
Nevertheless, these ideas have stimulated the first challeng- 
ing observations aimed at detecting the imprints of galaxy- 
IGM interplay. Adelberger et al. (2002, A02) obtained high 
resolution spectra of 8 bright QSOs at 3.1 < z < 4.1 
and spectroscopic redshifts for 431 Lyman-break galaxies 
(LBGs) at lower redshifts. By comparing the positions of the 
LBGs with the Lya absorption lines in QSO spectra, indeed 
they conclude that within ~ 0.5h~^ (comoving) Mpc of the 
galaxies little H i is present; on the contrary, between 1 and 
5h~^ Mpc an H i excess with respect to the IGM mean is de- 
tected. This simple interpretation might be at odd with the 
results of a VLT/UVES study of the Lya forest in the vicin- 
ity of the LBG MS1512-cB58 showing the opposite trend 
(Savaglio et al. 2002), i.e. an absorption excess close to the 
galaxy. 

Numerical simulations have also noticeable difficulties repro- 
ducing A02 results as discussed by Croft et al. (2002) and 
KoUmeier et al. (2002); however, these studies lack a self- 
consistent treatment of multiphase gas structure and/or out- 
flow dynamics. Here we revisit A02 results through SPH sim- 
ulations (Marri et al. 2003) that implement a new scheme for 
multiphase hydrodynamics and, more importantly, a phys- 
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ically meaningful outflow treatment; the full description of 
the code and of the tests made are given in Marri & White 
(2002). We then derive synthetic absorption line spectra 
along lines of sight randomly traced through the simulation 
box at 2 ~ 3 and compare them directly with A02 data to 
investigate galactic feedback effects on the IGM. 



2 SIMULATIONS AND DATA ANALYSIS 

We have performed hydrodynamic simulations for a ACDM 
cosmological model with f2o = 0.3, Q.a = 0.7, fit = 0.04 
and h — 0.7 km s~^ Mpc~^. The initial power spectrum 
is cluster-normalized (erg = 0.9); periodic boundary condi- 
tions are adopted. We first obtained a set of low-resolution 
runs (64'^ particles for both gas and dark matter) in a 7h~^ 
comoving Mpc cube; these runs serve as a guide for more 
computationally expensive runs and for testing purposes. 
We consider a first model where the IGM multiphase struc- 
ture and galaxy outflows are deliberately ignored and a sec- 
ond one which includes a description of both these physical 
effects. We will refer to these runs as sph and mfb , respec- 
tively, with the same meaning (and parameters) adopted in 
the description of the low-resolution ACDM test problem de- 
scribed in Marri & White (2002). The high-resolution run, 
on which the main results of the present analysis are based, 
is a 128^ particles simulation in a 10.5/i~^ comoving Mpc 
box. For this run we only studied the full mfb model. 
Softening lenghts are flxed both in physical and comoving 
coordinates as required in GADGET (Springel, Yoshida & 
White 2001). Gas (dark matter) gravitational softening is 
approximately &h'^ kpc (8/i~^ kpc) comoving and 3/i~^ kpc 
(4/i^^ kpc) physical in the high-res case and scales accord- 
ing to particle number and box size in the low-res run. In all 
runs we include the effects of a UV background produced 
by QSOs and flltered through the IGM, whose shape and 
amplitude are taken from Haardt & Madau (1996). 



2.1 Synthetic Lya forest spectra 

To allow a direct comparison of simulation results with ob- 
servational data we construct synthetic Lya spectra from 
simulation outputs at redshift z — 3.17 {z — 3.27) for the 
low-res (high-res) runs. In the simulation box we randomly 
trace 60 lines-of-sight (LOS) parallel to the x-axis, each of 
which is discretized into Npi^c = 1024 pixels. In order to 
assign to each pixel a value for the hydrogen density, tem- 
perature and peculiar velocity we perform a standard SPH 
smoothing using the 32 closest neighbour SPH particles to 
the LOS pixel position. The neutral hydrogen density is de- 
rived using the code CLOUDY94|^ on each pixel, adopting 
the same UVB as in the simulation. The transmitted flux 
due to Lya absorption in the IGM is oc e~^, where r is 
the optical depth along the considered LOS. The Hubble 
velocity vh varies in the range (0,«5J"^), where the max- 
imum value is set by the box size. The contribution to r 
at the observed frequency corresponding to vnik), where 
1 ^ ^ 1024 is the pixel index, is given by 



McDonald et al. (2000) 
mfb 643 
sph 64= 
mfb 1283 




Figure 1. Probability distribution function of transmitted flux at 
2 = 3.17 for the low-res simulations mfb (solid line) and sph (dot- 
ted) runs, and at z = 3.27 for the high resolution mfb run (dotted- 
dashed). The points represent the observational data (McDonald 
et al. 2000) at z = 3. 
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where Aa; is the pixel size, ctq is the Lya cross section, / 
is the Lya oscillator strength, Aq is the Lya wavelength, 
njj J (i), b{i) and v{i) = vnii) + Vpecii) are the H i number 
density, the Doppler parameter, and the total velocity in the 
pixel i, respectively; <l>g is a gaussian line profile and Vpec is 
the peculiar velocity. For additional discussion on this for- 
mula, see e.g. Ranch, Haehnelt & Steinmetz (1997a). 
To compare observations and simulations at best it is neces- 
sary to degrade the synthetic spectra to account for the un- 
certainties affecting the observed spectra. We perform such 
procedure through the following steps: i) continuum normal- 
ization at the highest flux value in each spectrum; ii) con- 
volution with the instrumental proflle; tii) sampling due to 
spectrograph spatial resolution; iv) addition of instrumental 
noise (see Ranch et al. 1997b; Theuns et al. 1998; McDonald 
et al. 2000; Petry et al. 2002). As most of the data we com- 
pare with are taken with the Keck/HIRES spectrograph we 
adopt the following instrumental characteristics: FWHM = 
8.0 km s~^, pixel spectral resolution AA = O.OSA and signal- 
to- noise ratio S/N=50. The low (high) resolution box has a 
wavelength extent of AAtox = 13.46A {AXtox ~ 20.83A); 
hence, degraded spectra in the two cases are made of 269 
and 417 pixels, respectively. 



3 RESULTS 
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As a general sanity check of the simulations we have flrst 
confronted the simulated statistical properties of the Lya 
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forest with the observed ones. The most obvious comparison 
involves the probability distribution function (PDF) of the 
transmitted flux. In Fig.^ we plot the PDF as a function of 
the flux (F) = for the low-res (both mf b and sph runs) 
and high-res simulations (mf b case only) and compare them 
with the observational data of McDonald et al. (2000). The 
general agreement is quite good through all the flux range 
and for all three cases. The discrepancy between simulations 
and data at high fluxes is probably due to the uncertainties 
in the data continuum fitting (McDonald et al. 2000; Croft 
et al. 2002). Surprisingly, it appears that the inclusion of 
multiphase and outflow physics, not considered in the pure 
sph run, does not affect the distribution in a sensible man- 
ner. In other words, galactic outflows leave the Lya forest 
unperturbed. This result is in agreement with that recently 
found by Theuns et al. (2002), who interpreted it as an indi- 
cation that galactic outflows tend to propagate preferentially 
in the voids leaving the Lya absorbing filaments virtually 
unaffected. As an additional check we have calculated the 
Doppler parameter, b, and the H i column density distribu- 
tion, and compared them with two QSOs observations at 
Zem = 3.38 and Zem = 3.17 by Hu et al. (1995). The val- 
ues of b and A'hi for each absorber in the synthetic spectra 
have been derived using the fitting programm AUTOVP. 
The experimental b distribution is well reproduced both by 
the sph and mfb (low-res + high-res) runs at similar qual- 
ity level (simulated and observed distributions peak both 
around fe ~ 20 km s~^). However, simulated spectra tend 
to slightly overproduce lines with 10 < 6 < 22 km s~^ and 
under-predict lines with 25 < & < 50 km s^^. A similar ef- 
fect has already been noted by Theuns et al. (1998). These 
authors propose a number of possible explanations for this 
behavior: physical (Hell reionization, radiative transfer ef- 
fects), numerical (resolution) and related to data analysis 
(fitting procedure). Our results seem to indicate that nu- 
merical artifacts should not be the dominant factor. As for 
the H I column density distribution, the mfb case seems to 
reproduce the observational data for log A^hi > 15.5 some- 
what better; for smaller values of A^hi the differences be- 
tween the two runs are negligible. 

In Fig.^ we show a temperature map from a slice through 
the high-res box a,t z = 3.27. Hot bubbles (T ~ lO'' K) of 
shocked gas produced by outfiows around the parent galax- 
ies are clearly apparent. Their sizes range from « 0.5 Mpc 
to « 2 Mpc , and their shape appears in some case 
rather jagged as a result of the interaction with the inhomo- 
geneous ambient medium. The volume filling factor of gas 
with temperature above 10^ K is 14%. The internal struc- 
ture of the bubbles can be inspected more quantitatively in 
Fig.^, where the comparison between sph and mfb is shown 
(low-res case) for various physical quantities (hydrogen den- 
sity, riH, ionization fraction, xhi , gas temperature, T, and 
peculiar velocity, Wpec) along the LOS through the center of 
the most massive galaxy M — 1.5 x lO^^Moh"^ in these 
simulations; the galaxy position corresponds to the density 
peak at a; « 6.8 Mpc . In both simulations the star for- 
mation rate in this galaxy is 5.6 A/0yr~^ (18 M©yr~^) for 
the mfb (sph ) models. A striking result of the comparison 
between the two models shows that, although outfiows are 
able to heat the halo/IGM gas up to high temperatures out 
to more than 1 Mpc from the galaxy, they do not seem 
able to modify its density structure in a sensible way. Hence 
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Figure 2. Temperature map in a slice through the high-res sim- 
ulation box {z = 3.27). The color bar shows the values of logT. 



the density in the surroundings remains high and close to 
that set up by the process of galaxy formation. Close to the 
galactic center the outfiow peak speed is about 130 km s~^, 
but this value rapidly decreases as kinetic energy is used to 
counteract the pressure of intergalactic accreting gas, rain- 
ing onto the galaxy at essentially the escape speed of the 
system, roughly 150 km s~^. The stalling radius is seen at 
the zero-crossing of Vpec, approximately 0.3 Mpc away 
from the outflow source. The relative insensitivity of the 
density to the SN energy injection can be interpreted as the 
fact that the outflow velocities are lower than the escape 
speeds and hence the flow is conflned by inflow. Also, we 
note that the outflow velocities we find are lower than those 
inferred by A02, 600 km s"^ from LBGs. 
Fig.^ shows the mean Lya forest flux averaged on pixels on 
different LOS at different distances, Ar, from the galaxy cen- 
ters in the simulation box, following the same experimental 
procedure as in A02: 



{F{Ar)) = 



1 



N{Ar) 



JV(Ar) 

■E 
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where N{Ar) is the number of pixels that fall in the bin 
Ar, and Fi is the transmitted flux in the pixel i. The syn- 
thetic data, normalized to the observational ones, are plot- 
ted only out to a distance Ar ~ 3.81 Mpc for the low- 
res simulations (Ar « 5.32 Mpc for high-res one) due 
to the periodic boundary conditions of the simulation. We 
flrst compare the two low-res sph and mfb cases. The dif- 
ference between these curves is negligible and on average 
the Lya absorption is similar: this is expected after the re- 
sults analyzed in Fig. 0. In addition both simulations are 
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Figure 3. Comparison between sph (dashed line) and mfb (solid 
line) low-res runs along the LOS through the center of the most 
massive galaxy (M = 1.5 X IO^^Mq/i"^) whose position corre- 
sponds to the density peak at a; ^ 6.8 Mpc h~^: a) hydrogen 
density, b) neutral hydrogen fraction, c) gas temperature, d) gas 
peculiar velocity. 



in agreement with the high-res one, thus ensuring that nu- 
merical convergence has been reached and resuhs are not 
affected by spurious effects. However, these models fail in 
reproducing the observed trend inside 1 Mpc . The sim- 
ulated Lya forest seems to be much more opaque than the 
observed one, with a maximum r discrepancy of about a 
factor 2.2. In order to assess if the disagreement could be 
reduced by a more efficient outflow clearing of the halos of 
smaller systems, we excluded the contribution of most mas- 
sive galaxies (M > 1.5 x IO^^/i^^A/q) to the flux. However, 
as it is seen from the Figure, this does not solve the prob- 
lem. As a flnal resort, we have randomized the position of 
galaxy centers in the box to account for the reported exper- 
imental error determination. To this aim we have added a 
gaussian random displacement to the redshift of each galaxy 
with r.m.s. cr^ = 0.002, which represents the 1-a error of the 
measure. This procedure improves the result, as now galaxy 
centers do not perfectly coincide with density peaks, but to 
an extent insufflcient to explain the data. 



4 CONCLUSIONS 

Motivated by the recent observational results of A02, we 
have studied, with the help of a set of cosmological simu- 
lations including star formation in multiphase gas and out- 
flows from galaxies, the effects of galaxy formation/activity 
on the properties of the surrounding Lya forest. Although 
on large scales our simulations can reproduce remarkably 
well a number of statistical properties of the IGM, they fail 
to predict the observed Lya flux increase in regions close to 



Figure 4. Mean Lya forest flux at different distances Ar from the 
galaxies in the simulation box. Synthetic results are shown for five 
cases: low-res sph (solid line) and mfb (dashed); low-res mfb with 
least massive (M < 1.5 X IO^^/i-^Mq) gal axies only (dotted); 
high-res mfb (dot-dashed); high-res mfb with randomized galaxy 
centers (long-dashed). The points are the observational results by 
A02 for LBGs. Normalization to the data at large Ar has been 
performed (see text). 



the galaxies themselves. The success can be ascertain to two 
concomitant effects: (i) outflows preferentially expand in re- 
gions of low-density (voids) thus preserving the filaments re- 
sponsible for the Lya absorbing network; (ii) the hot bubbles 
fill a relatively small fraction of the cosmic volume (« 14% 
in our simulations) . Support to the first hypothesis emerges 
also from an inspection of the velocity field in the surround- 
ing of the most massive galaxy in the simulation: the inflow 
of gas from low density regions is blocked by the bubble 
expansion, but it proceeds basically unimpeded along the 
fllaments (see Marri et al. 2003). 

Much more puzzling is instead the reason for the opacity 
excess (with respect to real data) we see in the simulation 
in the inner Mpc . Apparently, simulated outflows do 
not carry sufficient momentum to disperse the density peak 
created in the vicinity of the galaxy as a leftover of its forma- 
tion. Also, a large fraction of the outffow energy is used to 
counteract the infalling gas ram pressure which tends to pile 
up the gas into the galaxy. The temperature increase close to 
galaxies amplifies the magnitude of the coUisional ionization 
rate, which becomes larger than the equivalent photoioniza- 
tion rate for T 10^ K. In this case the H i neutral fraction 
xm is independent of gas density and remains at roughly 
the same level as in the general IGM. Therefore, coUisional 
ionization is not sufficient to balance the opacity increase 
induced by the density raise and the transmitted ffux drops 
accordingly to the latter. What are the possible alternative 
explanations for the observed ffux trend? The presence of a 
bias in the data produced by a preferential selection effect 
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of low Lya absorption LBGs has already been suggested 
by Croft et al. (2002). Another possibility is provided by 
photoionization, particularly if one recalls the recent results 
by Steidel, Pettini & Adelberger (2001), who detected flux 
beyond the Lyman limit (with significant residual flux at 
A < 912 A) in a composite spectrum of 29 LBGs at z = 3.4, 
a clue of a conspicuous escape probability of ionizing radia- 
tion from these objects. As mentioned above, the A02 data 
would require an optical depth (or equivalently a xm , as- 
suming a prescribed density profile) a factor ~ 2 — 3 smaller. 
If this can be achieved with the ionizing flux coming from 
galaxies must be proved with detailed radiative transfer cal- 
culations which are currently ongoing (Maselli et al. 2003) 
The first attempt using simplified analytical and/or post- 
processing techniques to account for this effect (Croft et 
al. 2002, KoUmeier et al. 2002) have yielded so far nega- 
tive answers, but a fully self-consistent, physically accurate 
description of the problem is awaited in order to draw a final 
conclusion. 
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